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Abstract - Granular matter exists out of thermal equilibrium, i.e. it is athermal. While conventional equi- 
librium statistical mechanics is not useful for characterizing granular materials, the idea of constructing a 
statistical ensemble analogous to its equilibrium counterpart to describe static granular matter was proposed 
by Edwards and Oakshott more than two decades ago. Recent years have seen several implementations of 
this idea. One of these is the stress ensemble, which is based on properties of the force moment tensor, and 
applies to frictional and frictionless grains. We demonstrate the full utility of this statistical framework in 
shear jammed (SJ) experimental states (T]|2), a special class of granular solids created by pure shear, which 
is a strictly non-equilbrium protocol for creating solids. We demonstrate that the stress ensemble provides an 
excellent quantitative description of fluctuations in experimental SJ states. We show that the stress fluctua- 
tions are controlled by a single tensorial quantity: the angoricity of the system, which is a direct analog of the 
thermodynamic temperature. SJ states exhibit significant correlations in local stresses and are thus inherently 
different from density-driven, isotropically jammed (JJ) states. 
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Introduction. A remarkable property of systems in thermal 
equilibrium is that the probability of occurrence of a micro- 
scopic state is known a priori through the universal Boltzmann 
distribution: P v = <?~P £v /Z(p). Here, (3 is the inverse tempera- 



C^ture and E v is the energy of microstate V. The Boltzmann dis 



tributions defines characteristics such as the relation between 
fluctuations and response [3|. By contrast, granular systems are 
intrinsically out of thermal equilibrium, and we lack a broad 
framework for describing their statistical properties [4|. The 
idea of constructing statistical ensembles to describe granular 
H systems originated in a proposal by Edwards |5 1 that the fluctu- 
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ations of slowly-driven, dense granular systems was controlled 
by the ensemble of "blocked" states, granular assemblies in 
static, mechanical equilibrium. The original Edwards ensemble 
uses free volume in a granular system as the analog of energy 
in a Gibbsian statistical framework. A more recent model is the 
force network ensemble (FNE) l6|[7), which has been useful in 
describing force fluctuations on fixed granular geometries. The 
stress ensemble is a generalization of the original Edwards idea, 
and is based on a "conservation" principle that arises from the 
constraint of local force and torque balance on every grain |8|. 



Merging of the original Edwards ideas with the stress ensemble 
approach leads to an analog of the Boltzmann distribution for 
granular solids (blocked or jammed states), where the role of 
temperature is played by two distinct quantities: (a) the com- 
pactivity which is conjugate to volume, and (b) angoricity, a 
tensor that is conjugate to the force-moment tensor of granular 
solids (9][T0). 

The applicability of the concept of angoricity as a 
temperature-like variable has been beautifully demonstrated 
in recent experiments using photoelastic disks JTTJ. In ear- 
lier work, the concept of angoricity was applied to analyze 
stress fluctuations in simulated assemblies of frictionless grains 
(9 12 1 . These granular solids, frictional assemblies in experi- 
ments, and frictionless in the simulations, shared the common 
feature that the solidification is density-driven. They, therefore 
fall within the rubric of the universal jamming phase diagram 



1 13 1, and can be viewed as zero-temperature limits of thermal, 



amorphous solids. 

A class of granular solids that are far-from-equilibrium, and 
not described within the original jamming scenario are shear- 
jammed (SJ) states. (Tl[2). These solids are created through 
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shearing without changing the density, which is a strictly non- 
equilibrium pathway for creating solids. They provide a unique 
opportunity for testing the stress ensemble with its full tensorial 
complexity in an ensemble that has no equilibrium analog. In 
this work we show that the stress fluctuations in SJ states are 
described quantitatively by the stress ensemble, which is the 
infinite compactivity limit of the generalized ensemble. This 
test of the stress ensemble is non-trivial since we demonstrate 
that stresses have non-trivial correlations in SJ states, unlike in 
density-driven granular solids. 

The paper is organized as follows: Section 2 provides a brief 
review of experimental SJ states, Section 3 introduces the ten- 
sorial formulation of the stress ensemble framework, Section 4 
presents results of testing the stress ensemble framework in SJ 
states, and finally, in Section 5, we discuss the correlations in 
SJ states, and its implications. 

Shear Jammed (SJ) states. Dry grains interact with purely 
repulsive, contact interactions. Without any cohesion between 
them, the only way to create a solid of dry grains is by ap- 
plying a load at the boundary. A common feature of most ex- 
perimental fTT|[T4) and numerical [15] techniques for creating 
jammed states is the application of isotropic pressure (we shall 
call IJ states from here on). The jamming transition is defined 
as the onset of mechanical stability, and usually associated with 
a particular packing fraction, (|)y, [13] that can be protocol de- 
pendent. IJ states can also be sheared to study their response 
[ 15 , 16 1. However, in these cases, shear stress should be consid- 
ered a perturbation to an already jammed state, and under large 
enough shear stress, IJ states can unjam into flowing states. 
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Fig. 1: Plot of the Pressure versus the packing fraction for all experi- 
mental states. Red: isotropic states with x ~ 0. Blue: states with x > 0. 
The experimental resolution for the pressure is Pq. below this thresh- 
old, there is uncertainty to weather a state actually has finite pressure, 
and is in fact jammed, or if the pressure is caused by friction with base 
plate. However, what is certain is that there are no isotropic states 
below (])/, for these states there is a sharp transition to finite pressure 
at §j. For states with shear stress, there is no sharp transition as a 
function of the packing fraction. 

In contrast, the observation of SJ states in experiments (TJ[2) 
and simulations JT7J has demonstrated that applying isotropic 



pressure is not a necessary condition for jamming. For a range 
of packing fractions below the isotropic (j) j for a given proto- 
col, jammed states can be created by applying shear only. The 
experimental SJ states, in two dimensions, are created through 
application of quasistatic, forward shear (TJ. Packing fraction 
does not play a significant role in determining the properties 
of SJ states. For example, it has been shown that pressure and 
shear stress do not depend monotonically on (j) as one would 
expect in IJ states, but rather monotonically scale with the frac- 
tion of force bearing grains in the system (TJ. The SJ states are 
qualitatively different from IJ states. Whereas IJ states can have 
non-vanishing pressure and zero shear stress, pressure arises 
as a consequence of shearing in SJ states |2[. SJ states have 
strongly anisotropic fabric and force networks and stress ten- 
sors. It has been shown that the SJ states emerge as a conse- 
quence of the percolation transition in the force network (TJ. 

Stress Ensemble. Over a decade ago, Edwards and 
Oakeshott [5| theorized that the dynamics of slowly driven 
granular materials is controlled by the statistics of mechani- 
cally stable configurations, know as blocked states. A simple 
thermal system is specified by the state variables E,V, and N. 
In constructing a statistical ensemble for athermal granular, the 
role of the Hamiltonian in a thermal system is replaced by the 
volume function of a granular packing. Known as the Edwards 
ensemble, all blocked states of the same volume V are assumed 
to be equiprobable (microcanonical hypothesis). An entropy is 
defined for a granular state with the state variables V and num- 
ber of grains N, S(V,N). The analog of the thermodynamic 
temperature can be defined by X(V) = dV /dS. X or the com- 
pactivity |5] is conjugate to the volume, and hence X = X(V) 
is also the primary equation of state. A microstate labelled by 
V, inside a granular packing with total volume V is found with 
a probability proportional to the canonical distribution 

P v = i«p[-Wv/X]. (1) 

Here, the volume function (analog of the energy E v s deter- 
mined by the set of grain positions w v — w({r,}). 

Recent interest has focused on conservation laws and the 
sampling of the phase space of blocked states, especially in the 
context of extending the ensemble framework to stiff but not in- 
finitely rigid grains. The basic aim of any ensemble approach is 
to predict the probability of occurrence of a microscopic state, 
given a set of macroscopic, measurable quantities such as the 
volume and the external stress. A static granular packing is 
ultimately characterized by the positions of the grains, {?,}, 
and the set of contact forces, {fij}. Specification of the forces 
is essential for defining static states of frictional grains, and it 
is therefore crucial to include the contact forces in a general- 
ized statistical description for infinitely rigid and stiff grains. 
Among the various formulations there are two different ways 
to incorporate {fij} in an ensemble, (1) the FNE [6 1 in which 
micro states are characterized only by {///}, and (2) the stress 
ensemble [8-10|[^| which considers {?,}, and {fij}. The latter, 
unlike the FNE, is not restricted to a fixed geometry 

'We emphasize that the stress ensemble described here is not on a fixed 
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In the stress ensemble, in addition to volume, granular 
blocked states are also characterized by the force-moment ten- 
sor, which is a function of both grain positions {r;;} and inter- 
granular forces {fij} 
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The constraint of local force balance leads to a conservation law 
for E ijsJITTJ, which states that the force-moment tensor cannot 
be changed by local intervention such as tapping. Changing E, 
while maintaining force balance on every grain, requires chang- 
ing the forces on a line of grains that spans the system [8 9 1 . For 
a large system, there is thus a conservation principle reminis- 
cent of energy conservation in non-dissipative systems. Taking 
tapping as an example for generating blocked states, the micro- 
canonical stress ensemble is defined by fixing the values of V 
and E. In addition to compactivity, the stress ensemble involves 
an intensive variable that is conjugate to E. This is the angoric- 
ity tensor, a -1 identified by Edwards and coworkers [5|. The 
relation between & and E is analogous to the relationship be- 
tween p and E in a thermodynamic system. 

In a canonical ensemble formulation J5J, the probability of 
obtaining a microstate v with volume w v and force moment ten- 
sor E v is given by a tensorial generalization of the Boltzmann 
distribution 
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Z(X,a) 
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exp (—a : E v ) , 



(3) 



where ":" means tensor contraction and the angoricity tensor 
0T 1 is defined by 



a k i 



In two dimensions, we have 



a : 
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In the above thermodynamic formulation, the entropy S(V,E) 
should be taken as the large M limit of the entropy of M grains. 
Since all blocked states are in force and torque balance, their 
force moment tensor E is always symmetric, hence a is also 
symmetric via eq. [4] It should also be pointed out that while 
both a and E transform like tensors under rotation, the argu- 
ment of the Boltzmann term in eq. [3] contains only terms that 
are rotationally invariant. 

Volume fluctuations in SJ states are less relevant than they 
are in IJ states. Analysis of SJ states shows that stress fluc- 
tuations and volume fluctuations are decoupled, and unlike IJ 
states, the packing fraction is not a state variable that deter- 
mines the jamming probability under the shearing protocol JT]. 
These observations suggest that the pure stress ensemble, ob- 
tained in the limit of X — > °° is the one relevant for describing 



geometry, but rather comprises of states with different contact and force net- 
works. 
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Fig. 2: (a) Schematic demonstration of a sub-region of size in and £ m 
inside a packing with (M, £at). Here we define subregions as any con- 
tiguous cluster of m-grains roughly circular in shape, (b) Eq.|6]predicts 
that when the distribution of local force moment tensor is multiplied 
by e a:J " (shown in inset for m = 15 and the 22— component of the lo- 
cal force moment tensor in SJ states), its logarithm should be related 
to the entropy of microstates S m , and for different a's this function 
should only differ by an additive constant. After a shift of each curve 
by Z22 calculated from the eq.[T2] we find the data collapse onto a uni- 
versal curve in the main figure. The functional form of the collapses 
data also shows good agreement with the postulated form of entropy 
(eq.|13}, given by the solid red line. 



the fluctuations of SJ states. For a subset m < M of the parti- 
cles, the pure stress ensemble |9j is defined by the canonical 
distribution 



P{%«) = exp [S m (E m )] exp (-a : t m ) . 



(6) 



The canonical distribution (eq. [6]l describes the stress fluctua- 
tion in a region containing the m grains that is in a "angoricity" 
bath due to the remainder of the system (see fig.[2ja)). S m {L m ) 
is the (microcanonical) entropy function which is defined by 
counting the number of microstates 



S m (t m ) oc log 



£ <*>v §(E„,-E V ) 



(7) 



In this formulation, we have allowed for the possibility of mi- 
crostates not being sampled equally by including a microscopic 
weight C0v for state v fl8) . The stress ensemble and the original 
Edwards ensemble can be extended to include this possibility 

(HID- 

Measuring Angoricity. In order to apply eq. [6] to the ex- 
perimentally generated SJ states, we need a measure of a. Un- 
fortunately, we have no 'thermometer' for directly measuring 
angoricity, and have to base its measurement on postulated or 
measured equations of state relating & to E. 

In IJ states, a scalar version of eq.[6]is applicable, and there is 
a scalar a that is conjugate to Y, the trace of the force-moment 
tensor. In jT2|, the stress ensemble was applied to analyze 
pressure fluctuations in simulated frictionless IJ states, and an 
equation of state was deduced from these measurements. More 
recently, an equation of state has been deduced from measure- 
ments in IJ states of frictional packings (TTJ. In both systems, 
a linear relation was found between the inverse angoricity and 
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Tm of an M grain system: a s» jf- . Since this linear relation 
worked quite well for pressure fluctuation in IJ states, we take 
a straightforward generalization of this as a postulate for the 
equation of state for SJ states. 



d <* M ( E, 
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Using this equation of state, each M-grain SJ state can be la- 
belled by an angoricity tensor, since the force-moment tensor 
is measured in the experiments. In the experiments, M ~ 1000, 
and one could wonder about this being large enough to be in 
the thermodynamic limit. We will address this question below 
through analysis of the convergence of the entropy function. 

Eq. [6] predicts the distribution of stress in a subregion con- 
taining m grains at the angoricity given by a. As has been 
pointed out before [19|, this Boltzmann-like distribution has 
a very special form. It has a term that depends purely on E, a 
term that depends purely on d, and the combination of these 
two variables appears only in the exponential. The stress en- 
semble would, therefore, predict that multiplying the distribu- 
tion by e wZ and taking its logarithm would yield a function 
that depends on E and distributions corresponding to different 
ds would differ only by an additive constant, or 



log 



P(t m )e &:± »> 



= S m (E m )-logZ(d) (9) 



We directly test these predictions as follows. All SJ states, 
regardless of the value of the shear strain or the packing frac- 
tion (j), are categorized by their global force-moment-tensor E^/. 
Through the equation of state (eqs.|8]l, the value of the angoric- 
ity tensor is determined. Then for each SJ state at a particular 
d, random contiguous clusters containing m grains (fig. |2ja)) 
are chosen to form an ensemble of subregions with different 
values of E m , giving a probability distribution (E m ) . To sim- 
plify the analysis, we avoid dealing with the multi-dimensional 
distribution function and instead, we analyze the reduced dis- 
tributions for each component of the local force moment tensor. 

Fig. |2jb) shows that the data collapse implied by eq.|9]works 
remarkably well for a subregion containing 15 grains. We find 
similar data collapse for the other components of the force mo- 
ment tensor, and for m as small as 4. For smaller m, there is 
more of a spread in the data, however, there is remarkably fast 
convergence as a function of m to a universal functional form, 
which we can interpret as the thermodynamic entropy S'(E^) 
for kl — 11,22, 12. The solid red line in fig. |2jb) is the func- 
tional form deduced from the assumed equation of state and the 
definition of of d (eq. |4| 



S m (±m) =AM log [det t M ] 



(10) 



Here, A = 1/2 is a constant of proportionality that is deter- 
mined from the best fit to the collapsed data. A few features of 
fig. |2|b) are of note: there is a small but systematic difference 
between the entropy deduced from the equation of state, and 
the form of the collapsed data, especially at small values of E^/. 
In spite of these small differences, the fact that we can collapse 



the data demonstrates that the Boltzmann-like distribution de- 
fined by an angoricity tensor (eq.[6]) works remarkably well for 
SJ states. 

The equation of state that we postulated is the analog of an 
ideal-gas equation of state relating temperature and energy. If 



the entropy implied by this equation of state (eq. 10 1 were to 
hold exactly for the SJ states, then it would imply that stress 
correlations in these states are shorter-ranged than the size of 
the m-grain subregions that we have analyzed. Below, we as- 
sume that the ideal-gas entropy holds, deduce an explicit scal- 
ing form for the distribution of stresses for m-grain subregions, 
and show that the ideal gas model fails to describe the distribu- 
tions in detail, and especially their scaling with m. To under- 
stand this observation, we analyze the variance of the compo- 
nents of the stress, and what it implies for stress correlations, 
and finally we discuss how to reconcile the success of the equa- 
tion of state with these observations. 

Since the SJ states are created with bi-axial shear the result- 
ing normal stress Eu m — E22,m is much larger (~ 10 times fTI) 
compared to the off-diagonal component E12 m- This justifies 
expanding the equation of state, eq.[8]to first order in 
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Then, reduced distributions of a particular force moment tensor 
component E,y can be obtained from eq. [6] and eq. [I0]by inte- 
grating out all other components mn ^ kl. For the 11 and 22 
components 

Pll,m(£ll,m) = / dl.22 dY, l2 P(E m ) 



a 



m/2+3/2 



pSredftll m) < , _0 'll E ll,nf 



P22,m(^22,n 



r(m/2 + 3/2) 
dL n dLn P(£ m ) 



(12) 



a 



m/2+3/2 
22 



r(m/2 + 3/2) 

Where we have introduced a reduced version of the entropy 



eq. 10 that depends on one force moment tensor component 

(£n, n ) (m/2+1/2) ; 



Sred{^22.m) = (^ll,m) W 2+1 / 2 ). 



(13) 



This form of the distribution for En,m and E22.H, at various 
d can be collapsed by defining the rescaled dimensionless vari- 
ables 



x = — anEn. m or 
m 

X = —<X22^22,m, 

m 



(14) 



yielding a gamma distribution 



_ m/2 + 3/2 /2+ i/2 —ikk 



(15) 
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The rescaled distributions from SJ states are compared to the 
prediction eq. 15 in fig. |3ja) for various local subregion sizes 
ranging from m = 5 to 71. This comparison shows that the 
theoretical prediction captures the mean and overall shape of 
the experimental distributions, but that theory and experiment 
differ significantly in the width of the distributions. 
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Fig. 3: (a) Probability distributions for the reduced variable defined in 
eq.|14| Solid red curves are theoretical predictions (eq. [13] with fitting 
parameter A = 1 /2) compared to SJ experimental data, m indicates 
number of grains in the subregion over which the local force moment 
tensor is defined, (b) Probability distributions for the reduced variable 
defined in eq. [16] Solid red curves are theoretical predictions com- 
pared to SJ experimental data, m indicates number of grains in the 
subregion over which the local force moment tensor is defined. 



There is no simple scaled form for the off-diagonal compo- 
nent Ei2,,„, but in the limit of eq. 1 1 or O42 — > 0, we can write 
a reduced distribution of Ei2. m in terms of rescaled local force 
moment components 



x = E 



12.1 



'22.1 



to obtain 



1^/2 + 3/2) 
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In fig.[3jb), we plot eq. [T7]for several m-values compared with 
SJ data. In this case, there is very good agreement between 
theory and experiment, including the variance. 

Discussion. Fig. [3|a) clearly shows that the variance of 
the local force moment tensor E is not well described by an 
ideal-gas type model,particularly for small m. To explore this 
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Fig. 4: log — log plot of the variances for each local force moment 
tensor component as function of m. The variances are scaled by arbi- 
trary stress scales to make them dimensionless. Dotted lines indicate 
the trend of m 1,5 . 



deviation further, we studied the scaling of the variance with 
the number of grains, m. As shown in fig. [4] the variances of all 
components of the stress scale as m 1 , with y m 1.5. We note that 
the variance is related to the integral of the two-point, spatial 
correlation function of the force-moment tensor of grains. The 
observed scaling, with an exponent larger than unity, can only 
arise if the correlations are long-ranged and do not decay within 
the region containing m grains. (20). Since our scaling is valid 
for m « 100, which corresponds approximately to a 10 x 10 
grain subregion, we can deduce that stress correlations do not 
decay away by 10 neighbors. The observed scaling of the vari- 
ance of stress components in the SJ states is in sharp contrast 
to the experimental observation in IJ states, where the variance 
scales as m fTT) . Interestingly, the variance of the scaled stress 
defined in eq. [16] scales as m indicating that this scaled stress 
does not have long-range correlations. For both IJ and SJ states, 
there is a well defined thermodynamic limit. SJ states converge 

with ygsEpy /(E) 

with v/((5EF)/(E} 

Since the probability of occurrence of SJ states is well- 
described by the Boltzmann-like distribution of eq. [6] we can 
use our experience with equilibrium statistical mechanics cal- 
culations to reconcile the adequacy of the ideal-gas like equa- 
tion of state with the presence of long-range correlations. It 
is known [20] that two-point correlations enter calculations of 
entropy (or free energy) through logarithmic corrections. It is, 
therefore, reasonable to expect that the entropy and the equa- 
tion of state are much less sensitive to the existence of cor- 
relations than the variance of the force-moment tensor. The 
long-ranged correlations in the local force moment tensor are 
indicative of long force chains that are clearly visible in SJ 
states |l][2T]. We are currently investigating the microscopic 
origin of the observed long-range correlations. 

In this work, we have demonstrated that the stress ensemble 
provides an excellent quantitative description of fluctuations in 
experimental SJ states. We show that the stress fluctuations are 
controlled by a single tensorial quantity-the angoricity of the 
system, which is a direct analog of the thermodynamic temper- 
ature. We show that the SJ states exhibit significant correlations 
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while IJ states converge faster 
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in local stresses in sharp contrast to IJ states. This observa- [21] Zhang J., Majmudar T. and Behringer R., Chaos, 18 
tion reinforces the conclusion [ 1 1 that SJ states are not merely (2008) 041107. 

density-driven jammed states created at a lower density. These 
states are inherently different from density-driven IJ states. 
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